AD-AlOt  797  ARMY  EN8INEER  WATERWAYS  EXPERIMENT  STATION  VICKSBUR9— ETC  P/9  8/11 

RAYLEX9H  WAVE  VELOCITY  MEASUREMENTS  USIN9  BROAD  BAND  FREQUENCY  — ETC(U) 
JUL  81  0  WILLIAMS 

UNCLASSIFIED  WES/MP/EL-91-3  NL 


( 


X 


Destroy  this  report  when  no  longer  needed.  Do  not  return 
it  to  the  originator. 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position  unless  so  designated 
by  other  authorized  documents. 


The  contents  of  this  report  are  not  to  be  used  for 
advertising,  publication,  or  promotional  purposes. 
Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of 
such  commercial  products. 


Unc  i  a. ;  n  i  fied 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Pete  Entered) 


[  REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

t.  REPORT  NUMBER 

Miscellaneous  Paper  Kb--'.--)  ^  ^ 

2.  GOVT  ACCESSION  NO 

3  RECIPIENT’S  CATALOG  NUMBER 

«•  TITLE  (And  Su6<«l.) 

RAYLEIGH  ^AVE  VELOCITY  MEASUREMENTS  USING 

BROAD  .BAND  FREQUENCY  SOURCES  - 

5.  TYPE  OF  REPORT  ft  PERIOD  COVERED 

-> 

1  Final  report 

6  PERFORMING  ORG.  REPORT  NUMBER 

7.  authorc;  . 

' . 1  {  >v  ■  , . 

Otis/Williams  j  j  ^  ,  J,.;  ■  //  /  _  >  \ 

6  CONTRACT  OR  GRANT  NUMBERf ■) 

9.  PERFORMING  ORGANIZATION  NAME  fHO  ADDRESS 

U.  S.  Army  Engineer  Waterways  Experiment  Station 
Environmental  Laboratory 

P.  0.  Box  631,  Vicksburg,  Miss.  39180  /  (, 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  ft  WORK JJN1T  NUMBERS 

Project  ^A762T30AT^2, 

Tasjjf  B/ETj .  ‘  ' 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Office,  Chief  of  Engineers,  U.  S.  Army  n 

Washington,  D.  C.  20311* 

12.  REPORT  OA-I*  , 

is.  NUMBER  OF  PAGES  ’  ■— *']T  / 

14.  MONITORING  AGENCY  NAME  A  AOORESSflf  dlfforwit  from  Controlling  Oltlce) 

15.  SECURITY  CLASS,  (of  thhVHfl&N)*”" 

Unclassified  | 

15a.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  thle  Report) 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  20,  If  different  from  Report) 


t0.  SUPPLEMENTARY  NOTES 

I  Available  from  National  Technical  Information  Service,  5285  Port  H< yal  , 

Springfield,  Va.  22151. 

19.  KEY  WORDS  (Continue  on  r.v.r..  etde  it  nlCMHfy  tend  Identity  by  block  number) 

Fourier  analysis.:  Wave  pr-'paratiiT. 

Rayleigh  waves 
Seismic  spnscrs 
Transducers 
Ye me  it;.  measurement 

20.  ABSTRACT  I ’Cootfnue  ea  r»rtrp#  tide  it  rrerreeeery  mod  identity  by  block  number ) 

Ray  i  e  i  rh  wave  y«l -~r  i  *.  i  t-.~  y.  ■  r<  *  measured  tv  r  br-ai  :  and  t.  ■:  c.e  exeit.fit.i  n  f 

nrnl  c.  sir  faces .  Tw.  n.  ise  :•••  urciv  were  used :  a  Libra' r  .is  in**  band-;  assed 

filtered  white  n>  ise  and  a  I. a.-.r.' ■  r  ir  i  {*•  ir. jeUr.es.  .'he  Rayleigh  wave  veb  ■- 
i’.ios  were  mputed  using  url*  :•  analysis  t  determine  phase  differ  “net's?  V 

*■  ween  r  ran.' iucer.'  !  i-.ateu  a  Ci  x<-  1  b  i  st  nrioe  apart.  t'.-nsurements  with  the  ir.pui- 

aj-'-e  act' -cl  with  the  v  i  *. -rat.  r  im.-«su rent-iil  s  and  showed  ■  .)  n  i  's  f  r 

v  :’.d  i :vif  ■!  > 


Unc lassi tied 

security  classification  ok  this 


DO  1473  EMTIO**  or  >  MOV  «*  IS  OBSOLETE 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  RAOEf»ll*n  Da it  EnianE) 


20.  ABSTRACT  (Continued). 

Rayleigh  wave  velocities  over  wide  frequency  bands.  The  impulsive  source  is 
recommended  as  a  more  field-expedient  method  for  computing  Rayleigh  wave  velocity 
as  long  as  the  transducers  used  do  not  saturate.  Best  results  can  be  obtained 
using  the  averages  of  repeated  impulses..'  — 


Unclassified 


PREFACE 


This  study  was  sponsored  by  the  U.  C .  Army  Corps  of  Engineers  in 
furtherance  of  DA  Project  No.  UAT62T30ATU2 ,  "Terrain/Operations  Simula¬ 
tions,"  Task  B/E3,  "Analytical  Techniques  for  the  Design  and  Application 
of  Sensors." 

This  study  was  conducted  by  the  IJ.  C.  Army  Engineer  Waterways 
Experiment  Station  (WES),  Vicksburg,  Miss.,  under  the  general  super¬ 
vision  of  Dr.  John  Harrison,  Chief,  Environmental  Laboratory  (EL), 

Mr.  B.  0.  Benn,  Chief,  Environmental  Systems  Division  (ESD),  and 
Mr.  J.  R.  Lundien,  Chief,  Battlefield  Environment  Team  (ESD).  The 
study  was  under  the  direct  supervision  of  Dr.  D.  H.  Cress.  Project 
engineer  was  MAJ  0.  Williams.  Personnel  making  significant  contribu¬ 
tions  were  SP5  G.  Radu,  ESD,  and  Mr.  M.  B.  Savage,  Instrumentation 
Services  Division.  SP5  Radu  assisted  in  analyzing  the  data  and  devel¬ 
oped  the  necessary  computer  programs.  Mr.  Savage  assisted  in  collection 
of  the  field  data.  This  report  was  prepared  by  MAJ  Williams. 

Commanders  and  Directors  of  WES  during  this  study  and  preparation 
of  the  report  were  COL  John  L.  Cannon,  CE,  and  COL  Nelson  P.  Conover, 

CE.  Technical  Director  was  Mr.  F.  R.  Brown. 

This  report  should  be  cited  as  follows: 

Williams,  0.  1981 .  "Rayleigh  Wave  Velocity  Measurements 

Using  Broad  Band  Frequency  Sources,"  Miscellaneous  Paper 
EL-81-3,  U.  S.  Army  Engineer  Waterways  Experiment  Station, 

CE,  Vicksburg,  Miss. 


CONTENTS 


Page 

PREFACE  .  1 

CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 


UNITS  OF  MEASUREMENT .  3 

PART  I:  INTRODUCTION .  ^ 

Background  .  ^ 

Purpose  .  6 

Approach  .  6 

PART  II:  CONCEPT .  7 

PART  III:  DATA  COLLECTION . 10 

Equipment . 10 

Test  Sites . 15 

Data  Collection  Procedures  .  15 

PART  IV:  DATA  ANALYSIS . 17 

Introduction  .  17 

Analysis  Results . 18 

PART  V:  CONCLUSIONS  AND  RECOMMENDATIONS . 30 

Conclusions . 30 

Recommendations  .  30 

REFERENCES . 31 

APPENDIX  A:  MEASUREMENT  OF  SPECTRAL  PROPERTIES  .  A1 

APPENDIX  B:  PHASE  RESPONSE  OF  L-1U  GEOPHONES  .  B1 

Shake  Table  Method . B1 

Impulse  Source  Method  .  B2 

Conclusions . B2 

APPENDIX  C:  SOIL-CEMENT  TEST  AREA . Cl 


2 


CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  use  1  in  this  report  can  be  con¬ 
verted  to  metric  (Sf)  units  as  follows: 


_ Multiply _ 

inches 

pounds  (force) 

pounds  (mass)  per  cubic  foot 


_ 

?  5.J» 
h.hhQ 
16.OI8I46 


_ To  Obtain _ 

millimetres 

newtons 

kilograms  per  cubic  metre 
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RAYLEIGH  WAVE  VELOCITY  MEASUREMENTS  USING 
BROAD  BAND  FREQUENCY  SOURCES 


PART  I:  INTRODUCTION 


Background 

1.  Rayleigh  wave  velocity  data  are  often  used  to  characterize  the 
near  surface  medium  rigidity  for  such  purposes  as  building  and  highway 
foundation  design  and  documenting  test  conditions  for  the  evaluation  of 
military  sensors  that  make  use  of  seismic  energy.  Because  of  the  close 
relationship  between  Rayleigh  wave  and  shear  wave  velocities,  the 
Rayleigh  wave  velocity  can  be  used  to  estimate  the  shear  modulus ,  and 

is  particularly  accurate  for  obtaining  such  data  near  the  surface. 

Within  several  metres  of  the  surface,  the  accuracy  of  the  Rayleigh  wave 
velocity  for  estimating  shear  modulus  is  somewhat  dependent  on  the  com¬ 
plexity  of  the  layering  of  the  medium  (Lundien  and  Nikodem  1973). 

2.  The  Rayleigh  wave  is  a  surface  wave  phenomenon  whose  degree 
of  influence  decreases  in  an  exponential  fashion  with  depth  beneath  the 
surface;  its  velocity  is  particularly  sensitive  to  medium  properties 
within  a  wavelength  of  the  surface.  In  order  to  characterize  the  shear 
modulus  very  near  to  the  surface  (i.e.  within  approximately  20  cm), 
which  is  often  required  in  studies  of  the  performance  of  battlefield 
seismic  sensors,  it  is  necessary  to  obtain  measurements  of  the  Rayleigh 
wave  velocity  for  frequencies  as  high  as  1000  Hz,  i.e.,  frequencies  con¬ 
siderably  higher  than  normally  associated  with  seismic  exploration  tech¬ 
niques  for  foundation  design.  In  situ  measurement  of  the  shear  modulus 
of  frozen  ground  is  an  example  of  a  situation  for  which  such  high 
frequencies  are  necessary  in  order  to  characterize  the  medium  within 
approximately  the  top  20  cm. 

3.  The  traditional  method  for  measuring  Rayleigh  wave  velocities  / 

has  been  to  use  a  vibrator  source-geophone  layout  where  the  vibrator 

source  has  a  driving  frequency  range  from  about  20  to  2000  Hz  (Heukelom 

U 


and  Foster  I960,  Chang  and  Ballard  1971,  Ballard  and  McLean  1975).  The 
procedure  for  measuring  the  Rayleigh  wave  velocity  has  consisted  of: 

a.  Estimating  the  Rayleigh  wave  velocity  for  the  medium  under 
investigation  for  the  lowest  frequency  available  from  the 
source  (usually  on  the  order  of  20  Hz)  and  spacing  the 
geophones  so  that  the  distance  between  geophones  is  on  the 
order  of  a  quarter  of  a  wavelength. 

(b.  Sinusoidally  exciting  the  medium  at  the  lowest  possible 
frequency  for  the  source  and  measuring  the  time  delay. 

The  frequency  is  then  increased  incrementally  (at  the 
discretion  of  the  operator)  and  the  time  delays  are  mea¬ 
sured  until  the  distance  between  geophones  is  on  the  order 
of  a  wavelength. 

c_.  Moving  the  geophones  closer  together  (approximately  a 

quarter  of  their  previous  spacing)  and  repeating  the  time 
delay  measurement  process.  The  velocities  at  each  fre¬ 
quency  are  then  computed  from  the  geophone  spacings  and 
the  measured  time  delays. 

1.  The  above  measurement  technique  is  subject  to  several 
limitations : 

a-  Field  data  collection  and  analysis  for  determining  veloc¬ 
ity  as  a  function  of  frequency  (with  a  frequency  resolu¬ 
tion  of  several  Hertz)  are  extremely  time-consuming  for  a 
broad  frequency  range. 

b.  If  the  operator  is  not  careful  to  properly  monitor  the 
wavelength  versus  geophone  spacing  relationships  identi¬ 
fied  in  3a  and  b,  large  errors  can  be  made  during  analysis 
because  of  failure  to  account  for  the  presence  of  more 
than  one  wavelength  between  geophones. 

c_.  At  high  frequencies  (exceeding  several  hundred  Hertz),  the 
wavelengths  are  so  short  that  procedures  3a  and  b  result 
in  very  close  spacing  of  the  geophones.  In  some  cases, 
the  geophone  dimensions  will  not  allow  the  geophones  to 
be  moved  sufficiently  close  together  to  ensure  that  less 
than  one  wavelength  occurs  between  geophones.  Errors  in 
locating  the  geophones  can  result  in  significant  errors 
in  velocity  calculation. 


5.  The  limitations  identified  above  in  measuring  Rayleigh  wave 
velocities  using  sinusoidal  waves  are  of  concern  in  site  characteriza¬ 
tion  for  military  seismic  sensor  studies  as  they  normally  must  be  car¬ 
ried  out  expeditiously  in  a  variety  of  terrain  conditions.  A  more 
appropriate  method  of  measuring  Rayleigh  wave  velocity  has  been 
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developed  that  makes  use  of  a  multi frequency  (broad  band)  source  and 
conventional  signal  processing  techniques;  this  report  presents  the 
results  of  applying  such  a  method  to  obtain  Rayleigh  wave  velocity  data. 

Purpose 

6.  The  purpose  of  the  study  reported  herein  was  to  demonstrate 
the  measurement  of  Rayleigh  wave  velocities  using  multi frequency  (broad 
band)  sources. 


Approach 

7.  The  source  of  vibrations  coupled  into  the  ground  was  program¬ 
med  to  provide  a  broad  range  of  frequencies.  Two  sites  having  distinctly 
different  constitutive  properties  were  selected  for  conduct  of  tests. 
Geophone  signatures  induced  by  the  broad  band  and  impulse  sources  were 
recorded  on  magnetic  tape,  digitized,  computer  processed,  and  velocities 
as  a  function  of  frequency  computed.  These  velocities  were  then  com¬ 
pared  with  those  measured  using  conventional  techniques. 


/ 
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PART  r I :  CONCEPT 


8.  Wien  the  surface  of  a  medium  is  disturbed  by  a  source,  wave 
motion  can  be  initiated  and  coupled  energy  propagated  outward.  The  mode 
of  propagation  can  be  Rayleigh  (i.e.  surface)  waves,  shear  (transverse) 
waves,  or  compressional  waves.  As  shown  by  Riehart,  Hall,  and  Woods 
(1970),  the  vibrational  component  induced  by  a  vertical  force  on  a  me¬ 
dium  surface  is  predominantly  Rayleigh  wave  motion,  particularly  at  dis¬ 
tances  on  the  order  of  several  wavelengths  or  more  from  the  source.  For 
the  Rayleigh  wave  motion,  the  amount  of  displacement  from  equilibrium  is 
denoted  by  £(r,t)  ,  where  r  is  the  range  from  the  source  to  the  mea¬ 
surement  point  and  t  is  the  time.  Since  the  propagation  is  assumed  to 
be  radially  outward,  the  amount  of  displacement  arising  from  a  single 
frequency  source  can  be  written  for  the  Rayleigh  wave  as: 

£(r,t)  =  A(r )  exp  ,1  ju  £t  -  yryyj  +  4  (l) 


where 

A(r)  =  amplitude  of  the  wave 

1* 

a)  =  angular  frequency 
v(w)  =  Rayleigh  wave  velocity 
0  =  arbitrary  phase  angle 

The  variables  in  parentheses,  r  and  w  ,  denote  a  functional  depen¬ 
dence  on  the  respective  variables.  The  quantity  sensed  by  the  geophones 
is  the  particle  velocity  £  ,  which  is  obtained  by  taking  the  time  de¬ 
rivative  of  £  .  The  resulting  expression  for  £  is  given  by: 


f 

r  t. 

1  r\  1*\ 

wA  exp  ,i  { 

3 

> 

1 

-P 

3 

J  +  9  +  2] 

If  two  geophones  are  located  at  distances  r^  and  r„  from  the  source, 
the  respective  particle  velocities  at  r  and  r0  are  given  by: 

.1  cl 


C(r1}t) 


=  wA  exp 


j 

_ I 

„  4 

(  OJ 

L  V  ( to )  J 

+  0  +  -> 

1 

2( 

(3) 
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;'(r0,t)  =  wA  exp  j 


where 


( 

(r  +  Ar) 

) 

\  0) 

_  1 

+  e  + 

3 

> 

1 

_> 

2i 

-  s 

r  —  r 

2  1 

(J»; 


(5) 


The  propagation  time  At  for  the  Rayleigh  wave  to  move  from  r^  to  r^ 
is  given  by: 


At  = 


Ar 

v(w) 


(6) 


As  may  be  seen  in  Equations  3  and  4,  the  above  propagation  time  results 
in  a  difference  in  phase  between  t! 
phase  difference  A9  is  given  by: 


in  a  difference  in  phase  between  the  wave  forms  at  r  and  r„  .  This 


AO  = 


to  A  r 
v(to) 

w  At  (to) 


(7) 


where 


A9 


=  phase  |  C(r1,t)|  -  phase  |c(r2,t)| 


(8) 


The  Rayleigh  wave  velocity  is  therefore  given  by: 


v(w)  = 


toAr 

A9 


(9) 


The  quantity  A9  is,  in  principle,  obtainable  from  the  measured  data. 
However,  mathematical  computation  of  A9  restricts  its  apparent  values 
to  a  domain  from  -it  to  +tr  (the  principal  values  of  A0  can  be  de¬ 
fined  to  be  0  <_  A9  <_  2tt  ;  however,  most  computer  algorithms  use  the 
principal  values  in  the  range  -it  <_  AO  <  it  ) .  Tf  the  propagation  time 
over  the  distance  Ar  exceeds  the  period  of  the  vibration  (but  is  less 
than  two  periods),  a  value  of  2tt  must  be  added  to  the  apparent  value 
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of  0  derived  from  mathematical  computation. 

9.  The  ambiguity  concerning  the  true  value  of  A0  discussed 
above  is  a  disguised  version  of  the  same  problem  encountered  in  the  con¬ 
ventional  procedure  for  Rayleigh  wave  velocity  measurement  in  which  the 
operator  spaces  the  geophones  no  more  than  one  wavelength  apart  (or  cor¬ 
rects  for  such  placement  in  subsequent  data  analysis).  As  shown  in  the 
next  paragraph,  the  apparent  ambiguity  can  be  removed  by  using  a  broad 
band  source  (rather  than  a  single  frequency  source)  for  Rayleigh  wave 
velocity  measurements. 

10.  For  a  broad  band  source  having  frequencies  ranging  from  the 
tens  of  Hertz  to  the  hundreds  of  Hertz,  the  phase  difference  at  each 
frequency  can  be  obtained  by  conventional  cross-spectrum  Fourier  trans¬ 
form  methods  (see  Appendix  A).  A  phase  difference  plot  as  a  function  of 
frequency  for  a  broad  band  source  might  appear  as  shown  in  Figure  1. 

For  frequencies  above  approximately  75  Hz,  an  integral  multiple  of  2tt 
must  be  added  to  the  apparent  phase  difference  if  relative  phase  accu¬ 
racy  is  to  be  preserved  over  the  entire  frequency  range  shown.  There¬ 
fore,  for  the  frequency  range  of  75  to  175  Hz,  the  total  phase  differ¬ 
ence  is  the  computed  phase  difference  plus  2tt  .  Similarly,  for  the 
frequency  region  of  175  to  275  Hz,  the  total  phase  difference  is  the 
computed  phase  difference  plus  Hit  .  The  velocity  at  each  frequency 
can  be  computed  using  the  total  phase  difference  and  applying 
Equation  9- 


Figure  1.  Phase  versus  frequency 


PART  J  7 1 :  DATA  COLLECTION 


11.  For  the  data  collection  effort,  two  kinds  of  sources  were 
used:  (a)  a  conventional  electromagnetic  vibrator  driven  b.  a  broad 

band  noise  source  (as  opposed  to  single  frequency  sine  waves)  and  (b)  an 
impulse  induced  by  a  hammer  blow.  The  data  collection  equipment  and 
procedures  (including  a  description  of  the  electromagnetic  vibrator)  are 
discussed  in  the  following  paragraphs.  Two  test  sites  were  selected  for 
the  data  collection  effort,  each  having  distinctly  different  media.  One 
test  site  had  a  deep  loess  soil;  the  other  had  relative  rigid  surface 
material  consisting  of  a  soil-cement  mixture  overlying  the  deep  loess 
soil.  The  test  sites  are  discussed  in  the  following  paragraphs  also. 


Equipment 


12.  The  equipment  configuration  for  the  Rayleigh  wave  velocity 
measurements  is  shown  in  Figure  2  for  the  vibrator  source.  Basic 
components  are : 


a. 

Pseudo-random  noise 

generator 

b. 

Band-pass  filters. 

c. 

Power  amplifier  for 

vibrator. 

d. 

Vibrator. 

e_. 

Geophones. 

f . 

Cascaded  amplifiers. 

£.* 

Tape  recorder. 

13-  A  description  of  these  components  follows: 

Pseudo-random  noise  generator.  A  General  Radio  Company 
random  noise  generator,  type  1390-B,  with  a  frequency 
range  of  5  Hz  to  5  MHz  was  used  as  a  broad  band,  uniform 
amplitude  (white  noise)  signal  source. 

b.  Band- pass  filters.  Krohn-hite  model  3202R  band-pass 
filters  were  used.  These  filters  have  two  channels  of 
input.  One  channel  is  operated  in  the  low- pass  mode  and 
the  other  channel  is  operated  in  the  high-pass  mode  to 
provide  the  desired  noise  bandwidth.  The  model  3202R 
will  function  as  a  band-pass  filter  with  an  attenuation 
rate  of  2h  db  per  octave  outside  the  pass  band.  The 
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Pseudo-random 
noise  generator 


Figure  2.  Equipment  configuration  for  Rayleigh  wave 
velocity  measurements  using  vibrator  source 
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Krohn-hite  filters  have  cutoff  frequencies  continuously 
adjustable  over  the  frequency  range  from  20  Hz  to  2  MHz. 

c_.  Power  amplifier  for  vibrator.  A  solid  state  Phase  Linear 
Model  1+00  direct  coupled  amplifier  was  used  to  power  the 
Goodman  vibrator.  The  amplifier  is  capable  of  producing 
700  W  of  signal  power  if  driven  to  full  output. 

<i.  Vibrator .  The  Goodman  vibrator  (50  lb  force*)  was  used 
as  an  energy  source.  The  vibrator  has  a  l6-ohm  coil  and 
has  a  power  rating  of  100  W. 

e_.  Geo  phones.  The  analog  signatures  resulting  from  the 

energy  source  were  measured  with  the  Mark  Products  model 
L-1U  geophones.  They  have  a  l+.5-Hz  natural  frequency  and 
are  commonly  used  in  seismic  refraction  work.  A  total  of 
four  L-1U  geophones  were  used.  The  amplitude  and  phase 
responses  of  the  L-1U  geophones  are  discussed  in 
Appendix  B  of  this  report. 

f_.  Cascaded  amplifiers.  Cascaded,  common-mode  rejection, 
operational  amplifiers  designed  and  constructed  at  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES) 
were  used  for  amplifying  the  geophone  signals.  These 
amplifiers  have  adjustable  gains  from  0.1  to  1000. 

£.  Tape  recorder.  An  PM  magnetic  tape  recorder,  Lockheed 
model  1+17,  which  records  and  reproduces  seven  channels 
of  data  at  1-7/8,  3-3/1+,  and  7-1/2  in. /sec,  giving  a  flat 
frequency  response  of  direct  current  (DC)  to  625  Hz,  DC 
to  1250  Hz,  and  DC  to  2500  Hz,  respectively,  was  used  for 
recording.  The  tape  recorder  also  records  in  the  direct 
mode.  During  these  tests,  all  signatures  were  recorded 
at  tape  speeds  of  7-1/2  in. /sec  and  band  limited  from  DC 
to  2500  Hz. 

lU.  A  photograph  of  the  recording,  monitoring,  and  amplifying 
instrumentation  is  presented  as  Figure  3.  The  L-1U  geophone  and  Goodman 
vibrator  are  shown  in  Figures  1+  and  5,  respectively.  Instrumentation 
identified  in  Figure  3,  but  not  discussed  as  basic  components  of  the 
data  collection  configuration,  are: 

a.  Tektronix  model  1+22  portable  oscilloscope. 

b.  Continental  Specialties  Corporation  electronic  counter 
(100  MHz  frequency  counter)  for  alignment  of  the  tape 
recorder . 


*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ment  to  metric  (Si)  is  presented  on  page  3. 
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Test  Site 


15-  The  general  test  area  was  located  at  WES  in  Vicksburg,  Miss., 
on  the  northeastern  most  part  of  the  installation.  The  in  situ  soils 
developed  from  loessial  deposits.  The  surface  soil  was  a  brown  to  dark 
brown  heavy  silt  loam  (CL)  and  the  subsoil  ranged  from  heavy  silt  loam 
(CL)  to  silty  clay  loam  (CL).  The  in  situ  medium  consisting  of  this 
loess  was  selected  for  one  of  the  two  test  sites  identified  in  para¬ 
graph  11.  Previous  measurements  of  shear  wave  velocity  of  the  loess 
medium  at  the  first  site  were  on  the  order  of  120  m/sec  within  several 
metres  of  the  surface.  The  measured  compression  wave  velocity  in  the 
upper  two  metres  was  about  350  m/sec.  The  second  site  (soil-cement 
site)  in  the  test  area  contained  a  mixture  of  loess,  sand,  and  cement 
placed  in  deep  loess  soil  in  a  pit  6.0  m  wide,  6.5  n  long,  and  0.92  m 
deep.  The  mixture  was  constructed  to  have  shear  and  compression  wave 
velocities  of  approximately  600  and  1000  m/sec,  respectively.  The  soil- 
cement  test  site  is  further  described  in  Appendix  C. 

Data  Collection  Procedures 

l6.  The  procedures  for  data  collection  consisted  of  the  following: 

a.  Equipment  setup.  Recording  instrumentation  and  source 
and  geophones  were  configured  for  data  collection. 
Geophones  were  placed  radially  outward  from  the  source 
(vibrator  or  impulse  derived  from  a  hammer  drop).  The 
relative  locations  of  the  source  and  geophones  are 
identified  in  Figure  6. 

b.  Source  activation.  For  the  vibrator,  the  driving  fre¬ 
quencies  were  band  passed  through  the  filters  identified 
previously.  The  frequency  bands  were:  25  to  100,  100 
to  200,  200  to  300,  300  to  UOO,  bOO  to  500,  500  to  TOO, 

TOO  to  1200,  1200  to  2^00,  and  25  to  2^00  Hz.  Once 
amplifier  gains  were  adjusted  appropriately  for  tape  re¬ 
cording,  30  sec  of  vibrator-induced  data  was  recorded. 
Immediately  following  the  vibrator  source  tests,  three 
successive  impulses  from  the  hammer  drop  were  recorded. 
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Layout  of  equipment  on  test  site 


PART  IV:  DATA  ANALYST V, 


Introduction 


17.  The  data  analysis  procedures  were  developed  to  emphasize  the 
following  points : 

a.  Demonstrate  the  measurement  of  the  Rayleigh  wave  velocity 
for  the  usual  frequency  region  for  such  measurements  (20 
to  approximately  1*00  Hz). 

b.  Evaluate  the  potential  for  extending  such  measurements  to 
frequencies  higher  than  1*00  Hz  and  identify  limitations 
in  such  an  extension. 

c_.  Define  alternative  instrumentation,  data  collection  pro¬ 
cedures,  and  analysis  techniques  that  might  be  employed 
to  improve  and  expedite  field  measurement  of  Rayleigh 
wave  velocities  using  broad  band  sources. 

18.  The  results  of  the  data  analysis  procedures  are  dependent  on 
the  methods  used  for  processing  the  geophone  signatures.  The  operations 
involved  in  such  an  analysis  consisted  of: 

a.  Digitizing  the  recorded  signatures. 

b.  Computing  the  Fourier  transform  and  cross-spectrum 
properties  of  the  signatures  from  any  two  geophone 
signatures . 

The  methods  used  for  cross-spectrum  analysis  to  compute  the  needed  phase 
difference  (Equation  7)  are  discussed  in  Appendix  A.  A  particularly 
useful  function  for  evaluation  of  the  accuracy  of  the  results  is  the 
squared  coherence  (Equation  A3,  Appendix  A).  The  squared  coherence  has 
values  ranging  from  0  to  1.  Values  of  squared  coherence  near  unity  are 
observed  when  the  phase  differences  between  two  signatures  in  a  narrow 
frequency  band  (the  band  width  depends  on  the  stability  of  the  signa¬ 
ture  characteristics  over  time)  are  consistent  for  repeated  samples 
selected  from  the  signature  record.  Therefore,  the  coherence  between 
the  geophone  signatures  was  computed  to  provide  at  least  a  qualitative 
evaluation  of  the  accuracy  of  the  estimate  of  the  phase  difference. 

19 •  The  results  reported  herein  were  obtained  using  the  following 
analysis  procedure: 
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a.  Analog  data  were  digitized  at  a  rate  of  5000  samples  per 
second. 

b_.  A  U096- point  Fourier  transform  was  computed  (for  a  0.8- 
sec  time  window). 

c.  Spectral  parameters  (Gx  ,  Gy  ,  Gxy  in  Appendix  A)  were 
averaged  over  four  successive  time  windows  (each  of 
0.8- sec  duration). 

Although  data  were  recorded  for  a  number  of  frequency  bands  (para¬ 
graph  l6b),  the  results  were  very  similar  for  successive  bands.  There¬ 
fore,  only  the  extreme  conditions  represented  by  the  lowest  band  (25  to 
100  Hz)  and  the  entire  band  (25  to  2^00  Hz)  are  presented  for  discussion. 


Analysis  Results 


Vibrator  source 

20.  Measured  phase  differences  between  the  signatures  for  geo¬ 
phones  1  and  2  for  the  loess  site  using  the  25-  to  100-Hz  random  noise 
signal  to  drive  the  vibrator  are  shown  in  Figure  7a.  The  phase  differ¬ 
ences  "flip"  from  +  tt  to  -  it  at  about  50  Hz  and  again  at  approxi¬ 
mately  135  Hz.  Although  the  driving  frequencies  were  filtered  outside 
the  25-  to  100-Hz  bank,  sufficient  energy  existed  in  the  driving  signals 
to  provide  continuous  phase  difference  information  to  well  above  100  Hz. 
The  squared  coherence  for  these  data  is  presented  in  Figure  7b.  When 
the  squared  coherence  is  continuously  greater  than  approximately  0.8 
over  a  bandwidth  exceeding  10  Hz,  it  can  generally  be  regarded  as  "good" 
in  the  sense  that  the  phase  difference  measurements  are  very  stable  from 
one  sampled  frequency  to  the  next.  The  coherence  is  quite  good  from 
about  20  Hz  to  approximately  300  Hz.  It  also  has  apparently  good  coher¬ 
ence  below  10  Hz.  However,  the  "good"  appearance  is  largely  an  artifact 
of  the  logarithmic  scale  for  which  the  first  cycle  (l  to  10  Hz)  occupies 
a  third  of  the  total  plot  so  that  the  high  coherence  levels  are  only 
several  Hertz  wide.  The  resulting  velocities  (derived  using  Equation  9 
for  which  the  total  phase  difference  is  used  by  including  the  factor  of 
2  it  increase  for  each  flip)  as  a  function  of  frequency  are  presented  in 
Figure  7c.  Figure  7c  shows  that  there  is  some  disperison  over  the 
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frequency  region  20  to  l60  Hz.  The  measured  velocity  is  about  120  m/sec 
at  20  Hz  fmd  drops  to  about  100  m/sec  between  100  and  l60  Hz.  These 
measurements  can  be  compared  with  those  tabulated  below,  which  were  ob¬ 
tained  on  the  same  soil  type  during  another  data  collection  effort  using 
the  single  frequency  (conventional)  Rayleigh  wave  velocity  measurement 
technique  discussed  in  paragraph  3. 


Rayleigh  Wave 

Soil  Medium  Frequency,  Hz  Velocity,  m/sec 


Soil  (loess)  100 

150 

200 

250 


125 

150 

125 

125 


The  latter  measurements  are  somewhat  higher  than  those  observed  using 
the  broad  band  noise  source  but  still  show  reasonable  agreement.  The 
Rayleigh  wave  velocities  were  obtained  on  somewhat  different  soil  con¬ 
ditions  (i.e.  the  soil  moisture  content  for  the  two  measurements  was 
not  identical)  and  the  differences  (i.e.  approximately  25  percent) 
are  considered  small  relative  to  their  range  for  naturally  occurring 
conditions  (i.e.  comparing  results  for  frozen  and  nonfrozen  condi¬ 
tions  for  soils  or  for  soil  and  rock). 

21.  The  spectrum  of  the  signal  induced  by  the  vibrator  driven  in 
the  frequency  region  25  to  100  Hz  is  presented  in  Figure  7d  and  Je  for 
geophones  1  and  2,  respectively.  The  amplitude  of  the  induced  vibra¬ 
tions  is  predominantly  in  the  region  70  to  about  150  Hz.  The  reasons 
for  the  low  amplitude  in  the  25-  to  50-Hz  region  are:  (a)  low  force 
levels  output  from  the  vibrator,  and  (b)  poorer  coupling  of  forces  to 
induce  vibrations  at  the  lower  frequencies  relative  to  the  coupling  at 
higher  frequencies. 

22.  Results  of  the  phase  difference  calculations  for  the  band 
25  to  21+00  Hz  are  presented  in  Figure  8a  for  geophones  2  and  3  for  in¬ 
stances  in  which  the  squared  coherence  exceeds  0.8.  The  squared  coher¬ 
ence  plot  for  this  band  is  presented  in  Figure  8b.  The  squared  coher¬ 
ence  is  generally  good  (exceeds  0.8)  across  the  range  20  to  1+50  Hz 
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except  in  the  frequency  region  30  to  U0  Hz.  The  resulting  velocities 
for  the  frequency  region  0  to  500  Hz  (for  which  the  squared  coherence 
exceeds  0.8)  are  presented  in  Figure  8c.  Although  the  data  are  not  pre¬ 
sented  herein,  for  frequencies  above  500  Hz,  the  squared  coherence  is 
so  poor  that  the  computed  velocities  are  highly  erratic  and  do  not  re¬ 
sult  in  a  continuous  curve  of  velocity  as  a  function  of  frequency.  This 
can  be  attributed  to  the  low  amplitude  of  the  signals  at  these  high  fre¬ 
quencies  and  to  the  erratic  phase  response  of  the  geophones  (see  Figure 
B1  in  Appendix  B).  The  low  amplitude  of  the  signals  at  the  high  fre¬ 
quency  (see  Figures  8d  and  8e)  is  primarily  due  to  the  attentuation  of 
the  propagating  waves  in  the  soil  (due  to  viscus  damping).  The  phase 
response  of  the  geophones  is  discussed  further  in  Appendix  B. 

23.  Results  of  velocity  computations  using  various  geophone  com¬ 
binations  (l  and  2,  2  and  3,  3  and  ^).on  the  loess  are  presented  in  Fig¬ 
ure  9a.  The  velocity  computations  were  made  using  the  25-  to  2h00-Hz 
broad  band  noise  source.  The  dependence  of  velocity  on  frequency  is 
about  the  same  for  each  pair  of  signals  used  in  the  computation.  How¬ 
ever,  there  is  a  20  to  30  percent  difference  in  the  calculated  veloc¬ 
ities  for  the  maximum  computed  velocity  (geophones  2  and  3  or  3  and  U) 
and  the  minimum  computed  velocity  (usually  geophones  1  and  2).  Differ¬ 
ences  in  velocity  derived  from  different  pairs  of  geophones  may  be 
attributable  to:  (a)  geometric  dependent  mixing  of  Rayleigh  wave  motion 
induced  by  the  source  with  other  motions  (such  as  background  noise  com¬ 
ing  from  a  direction  different  from  that  of  the  source,  reflections  of 
compression  and  shear  waves  off  of  deeper  layers,  etc.);  (b)  differences 
in  the  phase  response  of  the  geophones  used  for  recording;  and  (c)  non¬ 
homogeneities  in  the  media  between  the  pairs  of  geophones.  To  the  ex¬ 
tent  that  b  and  £  contribute  to  differences  among  geophone  pairs,  the 
most  representative  velocity  is  probably  best  determined  by  averaging 
the  three  velocity  curves.  (Similar  comparisons  of  measurements  from 
geophones  1,  2,  3,  and  U  were  made  for  the  soil-cement  site  and  the 
results  are  shown  on  Figure  9b;  also  see  comments  about  results  at 
frequencies  above  500  Hz  in  paragraph  25.) 

2U.  The  phase  differences  derived  from  geophones  1  and  2  on  the 
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Figure  9-  Effect  of  geophone  combination  for  measurements 
in  loess  and  soil  cement 


r~ 


soil-cement  site  for  the  25-  to  100-Hz  source-driving  frequencies  are 
presented  in  Figure  10a.  No  phase  flips  occur  over  the  frequency  range 
20  to  170  Hz.  As  shown  in  Figure  10b,  the  coherence  is  good  over  the 
frequency  range  tested  and  in  fact  good  to  290  Hz.  Since  the  distance 
between  geophones  1  and  2  is  the  same  for  the  soil-cement  site  as  for 
the  loess  site,  the  velocity  in  the  soil  cement  is  considerably  greater 
than  in  the  loess  as  evidenced  by  the  phase  flips  for  the  loess  (Fig¬ 
ure  7a).  The  resulting  measurements  of  the  Rayleigh  wave  velocity  over 
the  20-  to  170-Hz  frequency  range  are  presented  in  Figure  10c.  Again, 
the  low  amplitude  at  the  higher  frequencies  (Figure  lOd  and  lOe)  is  due 
primarily  to  filtering  out  the  driving  frequencies  to  the  vibrator 
above  170  Hz. 

25-  The  phase  differences  derived  from  the  25-  to  2l+00-Hz  source¬ 
driving  frequencies  are  presented  in  Figure  11a.  Phase  flips  are  ap¬ 
parent  at  approximately  200  and  700  Hz.  Both  the  phase  difference  plot 
and  the  squared  coherence  plot  (Figure  lib)  show  a  continuous  nature  to 
frequencies  as  high  as  1000  Hz.  The  presence  of  consistent  phase  dif¬ 
ference  information  to  such  high  frequencies  as  apparent  in  the  presen¬ 
tation  of  the  squared  coherence  (Figure  lib)  can  be  best  accounted  for 
by  the  fact  that  high  frequencies  are  attenuated  less  severely  with 
distance  in  the  rigid  soil-cement  medium  than  in  the  loess  medium. 
Therefore,  the  signal,  levels  of  the  propagating  waves  are  well  above 
the  recording  noise  level  at  the  higher  frequencies  for  the  soil-cement 
medium.  The  velocities  resulting  from  the  phase  difference  computation 
for  the  20-  to  2U00-Hz  band  are  presented  in  Figure  11c.  The  curve 
shows  the  velocity  increasing  as  the  frequency  increases  with  the 
velocity  at  high  frequencies  approaching  500  m/sec.  However,  as  shown 
in  Appendix  B,  the  difference  in  phase  responses  of  L-1U  geophones  be¬ 
comes  questionable  at  frequencies  over  500  Hz.  This  point  is  also 
demonstrated  in  Figure  9b,  as  the  velocity  plots  above  550  Hz  are  some¬ 
what  erratic.  An  alternate  transducer,  such  as  an  accelerometer,  would 
be  better  suited  for  making  velocity  measurements  at  such  frequencies. 
Nevertheless,  the  clear  interpretation  of  the  phase  flips  over  a  broad 
range  of  frequencies  (see  Figures  lid  and  lie)  demonstrates  the 
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d.  INDUCED  FREQUENCY  DOMAIN 
FOR  GEOPHONE  1 


Figure  11.  Results  for  soil 
cement,  25  to  2^00  Hz,  geo¬ 
phones  1  and  2 


*.  INDUCED  FREQUENCY  DOMAIN 
FOR  OEOPHONE  2 


advantages  of  the  broad  band  source  and  associated  signal  processing  for 
deriving  the  Rayleigh  wave  velocity,  particularly  at  high  frequencies. 
Measurements  taken  on  the  same  soil  type  (soil  cement)  using  the  conven¬ 
tional  single  frequency  Rayleigh  wave  velocity  measurement  technique 
are  tabulated  below: 

Rayleigh  Wave 

Soil  Medium  Frequency,  Hz  Velocity,  m/sec 


Soil  cement 

250 

360 

300 

375 

350 

375 

*400 

375 

Although  data  at  frequencies  above  UOO  Hz  for  the  conventional  single 
frequency  Rayleigh  wave  velocity  measurement  technique  are  not  presented 
here,  a  comparison  of  the  data  from  both  techniques  clearly  indicates 
that  the  broad  band  source  technique  has  potential  for  providing  more 
consistent  results  at  the  higher  frequencies. 

Impulse  source 

26.  In  field  implementation  of  measurement  of  Rayleigh  wave  veloc¬ 
ities  using  broad  band  sources,  use  of  an  impulse  source,  such  as  a  ham¬ 
mer  drop,  would  simplify  test  requirements  and  would  be  particularly  ex¬ 
pedient  for  military  application.  This  section  of  the  report  presents 
the  results  of  processing  the  geophone  signals  derived  from  hammer  drops 
to  obtain  the  Rayleigh  wave  velocity.  Figure  12a  presents  the  computed 
phase  differences  derived  from  a  single  hammer  drop  2  m  from  geophone  1 
for  the  loess  site.  The  phase  difference  curve  is  erratic  in  the  55-  to 
70-Hz  region  but  has  fairly  good  continuity  in  the  frequency  regions  20 
to  55  and  70  to  100  Hz.  Because  of  the  short  duration  of  the  impulse, 
it  is  difficult  to  get  an  accurate  measure  of  the  squared  coherence. 
However,  the  computed  squared  coherence  is  presented  in  Figure  12b  and 
shows  a  definite  drop  in  the  frequency  region  of  55  to  70  Hz. 

27.  The  Rayleigh  wave  velocity  computed  from  the  hammer  drop  is 
presented  in  Figure  12c.  In  the  55-  to  70-Hz  region,  the  velocity  devi¬ 
ates  dramatically  from  20  to  220  m/sec.  If  the  55-  to  70-Hz  frequency 
region  is  overlooked,  the  measured  Rayleigh  wave  velocity  derived  from 
the  impulse  agrees  reasonably  well  with  that  observed  for  the  loess 
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using  the  broad  band  vibrator  source  (compare  Figure  9a  with  Figure  12c). 
The  measured  velocities  for  the  impulse  source  are  somewhat  higher  in 
the  frequency  bands  20  to  55  Hz  than  for  the  vibrator  results  but  agree 
very  well  over  the  70-  to  l60-Hz  frequency  region.  The  lack  of  coher¬ 
ence  in  the  55-  to  70-Hz  frequency  region  may  be  related  to  the  dip  in 
the  amplitude  spectrum  in  this  region  (Figures  12d  and  12e).  Improve¬ 
ment  in  the  velocity  accuracy  over  the  20-  to  55-Hz  band  could  probably 
be  obtained  by  increasing  the  signal  amplitude  at  these  frequencies. 
However,  because  of  its  short  duration,  a  single  impulse  does  not  pro¬ 
vide  much  statistical  opportunity  for  averaging  cross-spectral  param¬ 
eters  (Appendix  A).  Impulses  may  provide  adequate  information  for 
deriving  the  Rayleigh  wave  velocity  curve  provided  repeated  impulses 
are  incorporated  into  an  averaging  procedure. 
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FART  V:  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


28.  Based  upon  the  study  reported  herein,  the  following  conclu- 
can  be  drawn: 

a.  The  limitations  in  measuring  Rayleigh  wave  velocities 
using  conventional  monofrequency  techniques  at  fre¬ 
quencies  higher  than  200  Hz  can  be  circumvented  by  using 
a  multi frequency  (broad  band)  source  and  special  signal 
processing  techniques  (paragraphs  7  and  19)-  Processing 
of  geophone  signatures  by  computation  of  the  Fourier 
transform  and  cross-spectrum  properties  (phase  differ¬ 
ence  and  squared  coherence)  will  indicate  where  velocity 
measurements  are  reliable  and  where  they  are  questionable. 

b.  The  results  of  the  impulse  source  method,  described 
herein,  compare  very  favorably  with  the  vibrator  source 
that  generates  frequencies  over  a  broad  region.  Impulses 
may  provide  adequate  information  for  deriving  the  Rayleigh 
wave  velocity  curve  provided  repeated  impulses  are  incor¬ 
porated  into  an  averaging  procedure  (paragraph  27). 

c_.  Rayleigh  wave  velocity  measurements  of  frequencies  above 
600  Hz  are  for  the  most  part  limited  by  the  measurement 
transducers  (paragraph  25).  The  accuracy  of  the  procedure 
described  herein  is  limited  only  by  the  instrumentation 
and  equipment  used  (especially  the  phase  response  of  the 
geophones ) . 


Recommendations 


29.  It  is  recommended  that: 

a.  Alternate  transducers  that  will  more  reliably  measure  high 
frequency  signals  (exceeding  several  hundred  Hertz)  be 
used  (i.e.  accelerometers)  in  place  of  geophones. 

b.  More  extensive  tests  be  conducted  to  improve  the  accuracy 
using  the  impulse  source  technique  where  repeated  impulses 
are  incorporated  into  an  averaging  procedure. 

c_.  Commercially  available  spectrum  analyzers,  capable  of  com¬ 
puting  the  Fourier  transform  and  spectral  properties  of 
multi  frequency  signals  within  seconds  after  the  data  have 
been  measured,  be  incorporated  and  used  along  with  the 
procedures  described  in  this  study  to  compute  phase  veloc¬ 
ity  for  Rayleigh  waves  propagating  in  terrain  materials . 
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APPENDIX  A:  MEASUREMENT  OF  SPECTRAL  PROPERTIES 

1.  This  appendix  discusses  the  measurement  of  phase  difference, 
wave  velocity,  and  the  squared  cotierence  of  signals  using  Fast  Fourier 
Transform  (FFT)  techniques. 

2.  The  geometric  relationship  addressed  in  this  discussion  is 
depleted  in  Figure  Al.  As  shown,  it  consists  of  sensors  spatially 
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Figure  Al.  Geometric  relation  between  sensors  and  incoming  plane  waves 

separated  by  a  soil  media  (soil/soil  cement).  The  energy  source 
(Goodman  vibrator)  produces  a  series  of  waves  that  are  assumed  to  be 
propagating  away  from  the  source  along  the  line  in  which  the  geophones 
are  placed.  For  the  purpose  of  this  study,  it  is  assumed  that  a  single 
mode  of  propagation  is  employed.  This  avoids  the  ambiguity  that  may 
occur  in  estimating  wave  velocities  when  several  vibrational  modes 
having  different  propagational  velocities  are  present,  hence,  leading 
to  different  velocity  measurements  for  each  mode  of  vibration  (Stabler, 
Baylot ,  and  Cress  1976). 

3.  When  energy  is  input  into  the  soil,  as  in  this  situation,  it 
is  spread  over  bands  encompassing  a  broad  range  of  frequencies.  When  a 
broad  region  of  frequencies  is  present  in  the  signatures,  a  plot  of 
phase  difference  versus  frequency  can  be  generated.  It  is  important  at 
this  point  to  recognize  that  the  speed  (velocity)  at  which  a  particular 
wave  travels  from  one  transducer  to  the  next  is  directly  related  to  the 
phase  difference  $  .  As  discussed  by  Stabler,  Baylot,  and  Cress  (1976) 
mathematically , 


Al 


<p  =  kd 


(Al) 


where 


k 

d 

f 

f/v 


then 


wave  number  =  2irf/v 

point  along  line  in  which  waved  are  propagating 

frequency 

velocity 

wavelength 


=  2Trfd/v 


To  eventually  find  the  velocity  v  of  the  wave,  the  quantity  is  given  by 


2vfd 


(A2) 


k.  The  process  of  computing  the  phase  differences  for  each  dis¬ 
crete  frequency  band  in  the  frequency  domain  has  been  made  easier  with 
the  advent  of  the  FFT.  It  is  now  practical  to  compute  the  whole  spec¬ 
trum  in  reasonable  computation  time.  Even  with  this  computation  tech¬ 
nique  (FFT)  of  the  phase  difference,  some  applications  of  this  approach 
may  produce  inconsistent  results  (Stabler,  Baylot,  and  Cress  1976).  For 
example,  the  calculated  phase  difference  between  two  signatures  can  be 
so  erratic  in  some  frequency  bands  that  it  is  impossible  to  identify 

which  estimates,  if  any,  are  reliable.  As  discussed  by  Hammon  and 

2 

Hannan  ( 197*0 »  the  coherence  parameter  y  r  can  be  directly  related  to 
the  accuracy  of  the  cross-spectral  estimate  (hence  phase  difference)  at 
a  particular  frequency.  Because  the  accuracy  of  the  phase  difference 
estimate  improves  rapidly  as  the  coherence  increases  toward  unity,  the 
coherence  parameter  can  be  used  to  determine  the  frequencies  that  can 
be  used  to  provide  the  best  estimate  of  phase  difference,  thereby,  pro¬ 
viding  improved  estimates  of  wave  velocities.  The  squared  coherence 
is  computed  using  the  results  of  the  FFT.  In  the  FFT,  the  spectral 
content  of  a  time  series  is  found  using  (Flohr  and  Cress  1979) 
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X,  =  At  Y  x  exp  ( —12 II  f'  ,  rnAt) 
k  in  k 


where 


X  =  complex  amplitude  associated  with  frequency  f 
K  K. 

Real  (X^)  =  amplitude  of  cosine  transform  at  frequency  f 
Imaginary  (X  )  =  amplitude  of  sine  transform  at  frequency  f 

K.  K 

At  =  time  between  samples 

N  =  total  number  of  samples 

x  =  amplitude  of  sampled  function  at  time  =  mAt 
m 


i  =  V-l 

f,  =  k/NAt 
k 

k  =  1...N 


Then  the  auto  spectrum  is  denoted  by 


,Jxk  Ak*Xk 


Where  the  astrisk  denotes  complex  conjugation.  The  cross  spectrum 
for  two  time  series  x(t)  and  y(t)  is  given  by 


Gxyk  ~  Xk*V 


Then  the  coherence  parameter  y  ^  is  given  by 


2  G  2 

xvk  = 


G  ,G  , 

xk  yk 

where  the  circumflex  (*)  indicates  an  average  value.  For  our  analysis, 
the  averaging  was  done  over  adjacent  frequency  values  using  a  centrally 
weighted  weighting  function.  For  example 


-  V  r 

M  ",  x(k+j-!i)W. 
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APPENDIX  1$:  PHASE  RESPONSE  OK  1,-1 U  OEOPHONES 

1.  The  purpose  of  this  appendix  is  to  present  data  from  tests 
performed  to  verify  the  accuracy  and  reliability  of  the  L-1U  geophones 
used  in  this  study.  The  verification  is  based  on  manufacturer  specifi¬ 
cation.  The  L-1U  geophones  are  produced  by  Mark  Products,  Houston,  Tex., 
and  have  a  natural  frequency  of  4.5  Hz.  They  are  commonly  used  in  seis¬ 
mic  refraction  work.  Figure  4  of  the  main  text  is  a  picture  of  an  L-1U 
geophone.  A  total  of  four  L-1U  geophones  were  used. 

2.  The  conventional  methods  of  verifying  geophone  reliability  or 
accuracy  are  by  (a)  using  a  vibration  pickup  shake  table  to  sweep 
through  the  desired  frequencies  at  a  constant  velocity  and  plotting  the 
phase  and  amplitude  of  a  single  geophone,  or  (b)  checking  the  phase  dif¬ 
ference  between  geophones  using  the  impuLse  source  technique  described 
in  paragraph  2 6  of  the  main  text.  These  methods  are  described  in  detail 
in  the  paragraphs  below. 

Shake  Table  Method 

3.  Each  L-1U  geophone  was  individually  placed  on  a  Mb  electronics 
Hamden,  Conn. ,  type  120,  vibration  pickup  shake  table.  This  particular 
table  was  designed  to  have  axial  resonance  above  2  KHz,  which  is  well 
above  the  working  range  of  the  pickup  to  eliminate  any  effect  on  the 
operation  of  the  geophone.  Each  geophone  was  tested  at  a  constant 
velocity  (0.1  in. /sec)  and  at  the  frequencies  and  sweep  rates  listed 
below: 

Frequency,  Hz  Sweep  Rate,  dcra-ie 'min 

2-10  1 
10-100  1-0.2* 

100-200  0 . 2 

200-1000  0 . | 

*  Sweep  rate  was  decreased  by  0.1  decade/ 
min  for  each  10-Hz  frequency  increase 

During  each  test,  a  tracking  filter  ban!  width  of  1  Hz  war  used. 
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b.  The  results  of  the  shake  table  tests  are  shown  in  Figures  B1 
and.  B2.  Figure  B1  (a  thru  d)  indicates  that  the  geophones  were  reason¬ 
ably  close  to  manufacturer  specifications.  For  each  plot  of  the  phase, 
the  curve  passes  through  90  deg  at  the  natural  frequency  (4.5  Hz)  of  the 
L-1U  geophone.  The  phase  plots  (Figure  Bl)  also  indicate  that  the  per¬ 
formance  of  the  geophones  was  "good"  up  to  approximately  195  Hz.  Each 
geophone,  except  geophone  1,  indicates  consistent  performance  for  ap¬ 
proximately  195  Hz  to  TOO  Hz.  Geophone  1  shows  consistency  from  ap¬ 
proximately  195  Hz  to  520  Hz.  The  amplitude  plots  for  each  geophone 
(Figure  B2  (a  thru  d))  show  an  initial  peak  of  about  0.60  V  at  approxi¬ 
mately  5  Hz ,  which  is  appropriate  for  an  underdamped  geophone  at  the 
stated  natural  frequency  of  the  geophones.  Other  erratic  peaks  occur  at 
approximately  the  same  high  frequency  ranges  listed  for  the  phase  plots 
in  Figure  Bl.  The  upper  limit  for  reliable  consistency  also  matches 
that  of  Figure  Bl.  This  test  also  indicates  that  the  data  analysis  and 
conclusions  in  the  main  text  were  reasonably  correct. 

Impulse  Source  Method 

5.  As  stated  earlier,  the  steps  of  this  test  are  described  in 
paragraph  26  of  the  main  text,  the  only  difference  being  that  the  sig¬ 
nals  from  the  hammer  drop  are  analyzed  with  respect  to  the  phase  dif¬ 
ference  of  two  geophones  rather  than  from  the  source  to  one  geophone. 

The  phase  difference  was  computed  and  plotted  (Figure  B3)  for  comparison 
of  geophones  1  and  2,  (Figure  B3a),  geophone  1  and  3,  (Figure  B3b) ,  and 
geophones  2  and  3  (Figure  B3c).  Data  from  geophone  b  were  incoherent  for 
this  test  and  were  not  included  in  this  analysis.  In  Figures  B3a  and 
B3b,  the  phase  difference  stays  near  0  rad  to  approximately  550  Hz  where 
a  significant  roll-off  is  indicated.  In  Figure  B3c,  the  phase  differ¬ 
ence  stays  near  0  rad  to  approximately  625  Hz  before  any  roll-off  occurs. 

Conclusions 


6.  As  a  result  of  these  test  analyses,  it  can  be  concluded  that 
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or  rbake  table 


Figure  B2.  Amplitude  plots  on  shake  table 
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APPENDIX  C:  SOIL-CEMENT  TEST  AREA 


1.  The  soil-cement  test  area  is  located  within  a  large  general 
test  area  at  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES), 
Vicksburg,  Miss.;  the  large  area  is  described  in  paragraph  15  of  the 
main  text.  The  soil-cement  area  is  a  soil  mixture  of  loess,  sarnd,  and 
cement  placed  in  a  pit  6.0  m  wide,  6.5  m  long,  and  0.92  m  deep.  The 
original  purpose  of  this  test  area  was  to  simulate  coral  or  frozen 
ground  conditions  by  fabricating  a  soil  mass  that  would  produce  shear 
and  compression  wave  velocities  of  600  and  1000  m/sec,  respectively. 

With  that  purpose  in  mind,  various  materials  were  gathered,  mixed,  and 
placed  in  test  cylinders  until  the  desired  mix  or  combination  of  mate¬ 
rials  was  determined.  The  final  product  that  was  machine-mixed  by 
volume  contained  2/3  loess,  1/3  masonry  sand,  6  percent  cement  (Portland 
Type  l),  and  10  percent  moisture.  The  soil-cement  mixture  was  placed  in 
the  pit  described  above  in  6-in.  lifts  and  compacted  using  commercial 
power  tampers.  Each  lift  was  compacted  to  an  average  density  of  125  pcf. 
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